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Initial indications for the production of a strongly coupled plasma 
in Pb+Pb collisions at ^s^v = 2.76 TeV 
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Results from first measurements of charged particle differential elliptic flow, obtained in Pb+Pb 
collisions at y/sJ^N = 2.76 TeV with the ALICE detector at CERN's Large Hadron Collider (LHC), 
are compared to those obtained for Au+Au collisions at ^snn = 0.2 TeV with the PHENIX detector 
at BNL's Relativistic Heavy Ion Collider (RHIC). The comparisons, made as a function of centrality 
(cent) or the number of participant pairs (JV part ) and particle transverse momentum pr, indicate 
an excellent agreement between the magnitude and trends for the flow coefficients V2(pr, cent). 
Analysis indicates that the averaged specific viscosity of the quark gluon plasma (QGP) produced 
in LHC collisions, is similar to that for the strongly coupled QGP produced in RHIC collisions. 
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a First results from Pb+Pb collisions at ^/snn = 2.76 

9 TeV, from CERN's Large Hadron Collider (LHC) [1, 2] 

10 have initiated the highly anticipated explorations of the 

11 the high temperature, high entropy density domain of the 

12 QCD phase diagram. At ~ 14 times the energy of RHIC 

13 collisions, these Pb+Pb collisions are expected to create a 

14 rapidly thermalized plasma of quarks and gluons (QGP) 

15 at temperatures higher than those currently accessible at 

16 RHIC. The reported hadron multiplicity in these Pb+Pb 

17 collisions is ~ 1584 (or 8.3 per participating nucleon pair 
ia -/Vpart) for the most central 5% of the hadronic cross sec- 

19 tion [1] - a factor of 2.2 increase over that observed in 

20 central Au+Au collisions at RHIC (t/sJw = 0.2 TeV). 

21 Thus, it appears that one now has a lever arm for probing 

22 the QGPs viscosity and other transport properties to de- 

23 tcrminc if they evolve from the strongly coupled plasma 

24 observed at RHIC [3-7], towards the more weakly inter- 

25 acting, gaseous plasma state expected at asymptotically 

26 high temperatures. 

27 In non-central heavy ion collisions, the spacial asym- 

28 metry of an initial "almond-shaped" collision-zone leads 

29 to flow. That is, partonic interactions in this collision- 

30 zone drive uneven pressure gradients in- and out of the 

31 reaction plane and hence, a momentum anisotropy of the 

32 particles emitted about this plane. At mid-rapidity, the 

33 magnitude of this flow is frequently characterized with 

34 the even-order Fourier coefficients; v n = (e'"' A ^) , n = 

35 2,4,..., where A(j) is the azimuth of an emitted hadron 

36 about the reaction plane, and brackets denote averaging 

37 over particles and events. 

38 Because they are known to be sensitive to various 

39 transport properties of the expanding hot medium [8- 

40 17], the differential Fourier coefficients V2 (N part ), V2(pt) 

41 and V2(N part ,pT) have been extensively studied as a 

42 function of collision centrality (cent) and hadron trans- 

43 verse momentum px, in Au+Au collisions at RHIC 

44 (v^w = 0.06 - 0.2 TeV) [18-26]. Indeed, consider- 

45 able effort is currently being devoted to the quantita- 



46 tive extraction of the specific shear viscosity r\j s (i.e. the 

47 ratio of shear viscosity rj to entropy density s) via com- 

48 parisons to viscous relativistic hydrodynamic simulations 

49 [16, 17, 27-34], transport model calculations [14, 15, 35] 
so and hybrid approaches which involve the parametrization 

51 of scaling deviations from ideal hydrodynamic behavior 

52 [7, 10, 13, 36, 37]. 

53 With the advent of detailed t>2(cent,pT) data for 

54 Pb+Pb collisions at the LHC (y/sNN = 2.76 TeV), an 

55 important question is whether these new flow data give 

56 an early indication for a significant difference in the vis- 

57 cosity of the QGP produced in RHIC and LHC collisions? 

58 Such a difference might be expected because, relative to 

59 Au+Au collisions at RHIC, the measured multiplicity for 

60 Pb+Pb collisions at ^snn = 2.76 TeV, suggests an ap- 

61 proximate 30% increase in the temperature of the QGP 

62 produced in LHC collisions. 

63 The influence of - on anisotropic flow is especially 

64 transparent in studies involving the flow coefficient scaled 

65 by the initial eccentricity of the collision zone V2 ( N v* Tt 'PT) 

&2 (-'''part ) 

66 as illustrated in Fig. 1. Here, results from hydrodynamic 

67 simulations (with the code of Dusling and Teaney [38]) 
ea are shown for two different viscosity values. For ^ = 0, 

69 Fig. 1 (a) indicates an essentially flat dependence for 

70 V2 ^(j^ rt '^ in li ne with the expected scale invariance of 

71 perfect fluid hydrodynamics. By contrast, Fig. 1 (b) 

72 shows that the introduction of a viscosity (- = 0.2) re- 

73 duces the magnitude of V2{N pax t,pr) and breaks the scale 

74 invariance of ideal hydrodynamics evidenced in Fig. 1 

75 (a). That is, there are substantial pT-dependent devi- 

76 ations away from the essentially flat Af part dependence 

78 observed in Fig. 1 (a). 

79 Figure 2 shows that these predicted scaling devia- 

80 tions are found in actual experimental data [37]. It 
si shows eccentricity-scaled values of V2a(j>t, N pa , T t) (ob- 

82 tained with factorized Kharzeev-Levin-Nardi [MC-KLN] 

83 model eccentricities [39, 40]) for several pr cuts. The 
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FIG. 1. (color online) Comparison of 1)2/22 vs. N paTt for 
several pr selections, obtained from perfect fluid (a) and vis- 
cous (b) hydrodynamic simulations of Au+Au collisions. For 
these calculation, a Glauber initial eccentricities are use in 
conjunction with a lattice-based equation of state [38]. 
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FIG. 2. (color online) Comparison of 1)2/22 vs. iV pa rt (a) and 
1)4/24 vs. iVpart (b) for several pr selections as indicated. The 
dashed curves indicate a simultaneous fit to the data in (a) 
and (b) [for each pr] [37]. The 1)2,4 data are from Ref. [26]. 



84 low-pr selections show small scaling deviations, i.e. they 
as are almost flat. However, the data points slope upward 

86 progressively (from low to high -/V par t) as the (pr) is in- 

87 creased, reflecting an increase in the scaling deviations 
sa with (pr). 

89 These eccentricity-scaling deviations reflect the effects 

90 of viscosity, as well as its attendant influence on the emis- 

91 sion distribution (/) on the freeze-out surface. This dis- 

92 tribution can be expressed as [9, 38]; 

dN y Q | 5 j _ ^ L | c ( PT \ 2 a \ (i) 

dyprdprdcj) ' y \ 2} / J ' 

93 where /o is the equilibrium distribution, Tf is the freeze- 

94 out temperature, C ~ s^fj arLC ^ a ^ s estimated to be 

95 [37]; r is the time scale of the expansion. Note that 

96 the factor 5f results [explicitly] from a finite shear vis- 

97 cosity and is known to dominate the calculated viscous 

98 corrections to V2(pr) for pr > 1 GeV/c due to its strong 

99 p\ dependence [38]. Thus, a significant increase in the 

100 value of 2 would not only serve to decrease the magnitude 

101 of " 2 e ^jv rt 'f) T " > but would also magnify the eccentricity- 

103 scaling deviations, especially for pr > 1 GeV/c. 

104 Figures 1 and 2 show that a simple way to test for 

105 a change in ^ for two different data sets, is to compare 

106 their respective eccentricity-scaled anisotropy coefficients 

107 HaljWiPr) d v a (n t , P T) t if th diffcr _ That - g 

108 a significant -j difference would not only lead to different 

109 magnitudes, but also to very different p^-dependent cur- 
uo vaturcs for the eccentricity-scaled coefficients from each 



111 data set. If the -/Vp art dependence of £2,4 is the same for 

112 both data sets, then the test can be made more simple 

113 by directly comparing the flow coefficients ^(cent, pr)- 

114 Indeed, the calculated MC-KLN initial eccentricities for 

115 the two reactions are very similar as shown in Fig. 3 (b). 
lie The same trend is observed for Glauber initial eccentric- 
117 ities which are smaller than the MC-KLN values. The 
us ratios in Fig. 3 (b) are a little larger than unity due to 

119 the larger size of the Pb nucleus. However, for the same 

120 centrality, they are w 1 as also noted in Ref. [2]. 

121 The flow results recently reported in Ref. [2] have also 

122 indicated a strong similarity between the elliptic flow co- 

123 efficients D2(cent,pr) obtained by the ALICE collabora- 

124 tion for Pb+Pb collisions at ^/s^n = 2.76 TeV and those 

125 obtained by the STAR collaboration for Au+Au colli- 
de sions at ^/snn = 0.2 TeV. Given that the differences be- 

127 tween the Glauber-based initial eccentricities for Au+Au 

128 and Pb+Pb collisions are small for the same centrality 

129 selection (cf. Fig. 3 and Ref. [2]), the measured flow 

130 coefficients for both data sets can be directly compared 

131 to test for a viscosity difference. 

132 A comparison of V2(pr) for several centrality selections 

133 from the PHENIX [26] and ALICE [2] data sets, is shown 

134 in Fig. 3 (a). The comparison shows good agreement 

135 between the magnitudes and trends for both data sets, 

136 indicating a strong similarity between the viscous cor- 

137 rections to v-^ipr) in Pb+Pb (^/sjvjv = 2.76 TeV) and 

138 Au+Au (y/sNN = 0.2 TeV) collisions [41]. Parenthet- 

139 icaily, an exact agreement between the magnitudes of 
ho both data sets is not to be expected because the ALICE 
mi measurements were obtained via the 4-particlc cumulant 
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FIG. 3. (color online) Comparison of V2 vs. pr for several cen- 
trality selections as indicated (a). The ALICE and PHENIX 
data are from Refs. [2] and [26] respectively. The ratio of the 
initial eccentricity for Pb+Pb and Au+Au collisions is shown 
as a function of A/p art in panel (b). 



«2 method [42] while the PHENIX measurements were ob- 

143 tained via the event plane method, albeit with a sizable 

144 Ary-separation between the event plane and the detected 

145 hadrons [20]. These different measuring techniques re- 

146 fleet different associated eccentricity fluctuations which 

147 manifest as a small difference in the magnitudes of the 
us two data sets. This difference is illustrated in Fig. 4 

149 where we show the ratio of the PHENIX «2{2} measure- 

150 ments to STAR's four particle ^2{4} measurements. The 

151 ratios show the expected 9-12% difference (csentially in- 

152 dependent of pt) due to the larger inherent fluctations 

153 for the v 2 {2} measurements [43, 44]. This difference does 

154 not alter the arguments nor the conclusions which follow. 

155 The observed agreement between the ^(pt) data from 

156 both the LHC and RHIC implies that the observed in- 

157 crease of the pr-integrated v-i (from RHIC to the LHC) 

158 [2], can be simply explained by an increase in the (pr)- 

wo As in Refs. [7, 37] the deviations from eccentricity- 

161 scaling have been used to characterize the magnitude 

162 of the viscous corrections to v ^ n ^^t) an( j ^(AWt) 

£2 (.-"part J E2(-lV pa rt) 

163 [10, 13, 36, 45] by a Knudsen number (K = X/R) 

164 parametrization, where A is the mean free path and R is 

165 the transverse size of the system obtained from the same 

166 Glauber-based calculations used to determine S2(N palt )- 
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FIG. 4. (color online) Comparison of PHENIX's v 2 {2} vs. p T 
and STAR's «2{4} vs. pr for several centrality selections as 
indicated. The STAR and PHENIX data are from Refs. [2] 
and [26] respectively. The dotted and dashed lines indicate 
ratios of 1.0 and 1.15 respectively. 



In turn, the extracted Knudsen number provides an esti- 
mate for the specific viscosity of the QGP; 



1 

•s 



XTc s = (RKTc s ), 



(2) 



169 where c s is the sound speed estimated from lattice calcu- 

170 lations [46] for the mean temperature T . The agreement 

171 between the LHC and RHIC data shown in Fig. 3 (a) 

172 and in Fig. 2 of Ref. [2], indicate very similar viscous 

173 corrections and thus, a similar - range for the plasma 

174 produced at higher temperatures in Pb+Pb collisions at 

175 y/s^N = 2.76 TeV. In Ref. [37] the estimate 4/T§ - 1 - 2 

176 was obtained for the K values extracted using MC-KLN 

177 and MC-Glauber eccentricities [respectively] in central 

178 and mid-central Au+Au collisions (^/sjvlv = 0.2 TeV) 

179 for the mean temperature T = 220 ± 20 McV [47]. 

180 The similarity between the ^ values for the plasma pro- 
lsi duced in RHIC and LHC collisions can be understood in 

182 the framework of Eq. 2, via the following simple estimate 

183 for the Knudsen number [48, 49]; 



K 



J_ 
RT 



(3) 



where the magnitude of (3 depends primarily on whether 
the plasma is strongly or weakly coupled (for a weakly 
couple plasma, /3 ~ 36/8.144<7 4 ). Substitution of the es- 
timate for K into Eq. 2 shows that very little change in 
j would result if the coupling strength of the plasma re- 
mains essentially the same for two different mean temper- 
atures, i.e. the mean sound speed does not show a strong 
temperature dependence over the range of interest. Note 
that a similar argument applies for the comparison of 
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RHIC differential v-i data over the beam collision energy 
range ^snn = 0.062 — 0.2 TeV, where -^(pt, cent) has 
been observed to be approximately constant for Au+Au 
collisions [21]. Here, an important difference is that the 
associated temperature change is relatively small. 

In summary, we have made detailed comparisons be- 
tween measurements of charged particle differential el- 
liptic flow obtained in Pb+Pb collisions at ^/snn = 
2.76 TeV, and those obtained for Au+Au collisions at 
^/sWn = 0.2 TeV with the PHENIX detector at RHIC. 
The comparisons indicate an excellent agreement be- 
tween the magnitude and trends for the flow coefficients 
V2{pti cent). Our analysis indicates that the averaged 
specific viscosity of the QGP produced in LHC collisions 
is similar to that for the strongly coupled QGP produced 
in RHIC collisions. Therefore, a strong indication for an 
evolution toward a more weakly interacting plasma has 
not been exhibited. It will be most interesting to investi- 
gate whether or not this conclusion is further supported 
by detailed viscous hydrodynamic calculations, as well as 
more detailed differential flow measurements at the LHC. 
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